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Abstract 

Presented  is  a  quantum  lattice  gas  for  Navier-Stokes  fluid  dynamics 
simulation.  The  quantum  lattice-gas  transport  equation  at  the  micro¬ 
scopic  scale  is  presented  as  a  generalization  of  the  classical  lattice-gas 
transport  equation.  A  special  type  of  quantum  computer  network  is  pro¬ 
posed  that  is  suitable  for  implementing  the  quantum  lattice  gas.  The 
quantum  computer  network  undergoes  a  partial  collapse  of  the  wavefunc- 
tion  at  every  time  step  of  the  dynamical  evolution.  Each  quantum  com¬ 
puter  in  the  network  comprises  only  a  few  qubits,  which  are  entangled  for 
only  a  short  time  period.  A  Chapman-Enskog  type  analysis  of  the  quan¬ 
tum  computer  network  indicates  that  the  total  system  of  qubits  behaves 
exactly  like  a  viscous  lattice-gas  fluid  at  the  macroscopic  scale.  Because 
of  the  quantum  mechanical  nature  of  the  scattering  process,  superposition 
of  outgoing  collisional  possibilities  occurs.  The  quantum  lattice  gas  obeys 
detail  balance  in  its  collisions  and  is  therefore  an  unconditionally  stable 
algorithm  for  fluid  dynamics  simulation. 


1  Introduction 

Prior  to  the  advent  of  digital  computing  in  the  late  1940’s,  analog  computers 
held  much  promise.  An  electrical  circuit  can  be  constructed  to  simulate,  say,  an 
underdamped  oscillator  governed  by  a  second  order  differential  equation.  For 
example,  an  electrical  circuit  can  drive  a  trace  on  an  oscilloscope  mimicking  the 
vertical  motion  response  of  a  fast  moving  automobile  with  poor  shock  absorbers 
after  passing  over  a  speed  bump — one  continuous  physical  system  configured  to 
behave  just  like  another  continuous  physical  system.  Today,  after  five  decades 
of  digital  computing,  history  may  repeat  itself  in  the  sense  that  it  may  again  be 
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worthwhile  to  build  “analog”  computers — for  example,  a  quantum  mechanical 
spin  system  configured  to  behave  just  like  a  Navier-Stokes  fluid. 

The  purpose  of  this  paper  is  to  show  how  to  do  this.  We  show  how  to 
arrange  a  network  of  small  quantum  computers  so  that,  taken  as  a  system,  the 
qubits  within  the  network  mimic  the  behavior  of  a  system  of  massive  quantum 
particles  moving  and  colliding  on  a  discrete  spacetime  lattice.  This  discrete 
quantum  particle  system  is  termed  a  quantum  lattice  gas  and  the  associated 
quantum  computer  network  is  called  a  lattice- gas  quantum  computer. 

Over  a  decade  ago,  classical  lattice  gases  were  found  that  behave  like  a  vis¬ 
cous  Navier-Stokes  fluid  at  the  macroscopic  scale  [1,  2].  In  this  paper  we  show 
that  a  quantum  lattice  gas  does  too.  The  prediction  of  the  quantum  lattice  gas’ 
macroscopic  equations  of  motion  is  achieved  by  a  generalized  Chapman-Enskog 
analysis.  A  property  of  the  quantum  lattice  gas  (when  used  as  a  numerical  al¬ 
gorithm  to  implement  a  probabilistic  lattice  gas  on  a  classical  computer)  is  that 
continuous  macroscopic  fields  for  the  mass  and  momentum  densities  are  directly 
obtained.  There  is  no  need  for  either  ensemble  averaging  or  coarse-grain  space- 
time  averaging,  which  are  computationally  expensive  in  a  classical  lattice-gas 
simulation.  This  computationally  useful  property  of  the  quantum  lattice  gas 
arises  because  it  models  the  discrete  particle  system  directly  at  the  mesoscopic 
scale,  avoiding  noisy  fluctuations  while  retaining  detailed  balance  in  the  local 
particle  collisions  [3].  Detailed  balance  is  satisfied  because  of  the  unitary  action 
of  the  collision  operator  as  it  causes  quantum  mechanical  superpositions  of  out¬ 
going  particle  configurations  at  each  site  of  the  spatial  lattice.  Consequently, 
the  quantum  lattice  gas  is  unconditionally  stable  as  a  numerical  algorithm. 

We  calculate  the  single-particle  distribution  function  analytically  for  a  quan¬ 
tum  lattice-gas  system  at  local  equilibrium.  The  analytical  prediction  is  that 
it  has  the  same  form  as  the  single-particle  distribution  function  of  a  classical 
lattice-gas  system,  which  also  obeys  the  principle  of  detailed  balance.  We  ver¬ 
ify  this  prediction  through  numerical  simulation  of  a  two-dimensional  quantum 
lattice  gas,  which  is  a  straightforward  generalization  of  the  classical  FHP  lat¬ 
tice  gas  [2].  For  comparison  purposes,  we  also  include  results  from  a  classical 
FHP  simulation.  In  the  low  Mach  number  incompressible  fluid  regime,  there  is 
excellent  agreement  between  the  theoretical  prediction  and  the  numerical  data 
for  the  single-particle  occupation  probability. 

2  Review 

There  are  new  possibilities  and  limitations  that  arise  in  computing  if  we  use 
the  principle  of  quantum  mechanical  superposition  of  states  [4,  5,  6,  7,  8,  9]. 
In  quantum  computing  a  two-level  quantum  bit  represents  the  smallest  unit 
of  information  which  may  be  in  a  superposition  of  the  discrete  states  |0)  and 
ll}.1  An  example  of  the  physical  embodiment  of  a  qubit  is  the  z-component  of 

1 A  qubit,  «)  =  ajO)  +  /3|1),  has  an  amplitude,  a,  of  it  being  in  the  zero  state,  |0),  and 
another  amplitude,  ft,  of  finding  it  in  the  one  state,  |1).  The  probabilities  add  to  unity: 
(0|0)  -(-  (1 1 1)  =  (<i\q)  so  the  complex  coefficients  are  constrained  by  |a|2  +  |/3|2  =  1. 
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a  nuclear  spin  in  an  atom  in  a  uniform  external  magnetic  field.2 

An  open  issue  for  quantum  computing  is  whether  entangled  qubit  states  (of 
many  qubits,  much  more  than  two)  can  be  isolated  from  the  surrounding  envi¬ 
ronment  for  delicate  quantum  algorithmic  steps  to  be  completed.  Using  quan¬ 
tum  mechanical  superposition  among  qubit  states  to  speedup  a  computation  by 
simultaneously  encoding  many  possibilities,  an  approach  termed  quantum  paral¬ 
lelism,  is  generally  considered  the  primary  virtue  of  quantum  computation.3  Yet 
uncontrolled  coupling  with  the  surrounding  environment  causes  decoherence  of 
the  qubit  states  and  the  virtue  is  lost — quantum  parallelism  levies  a  high  price 
for  coherence  of  the  quantum  computer’s  wavefunction.  This  has  spurred  the 
development  of  scalable  quantum  error  correction  techniques,  considered  cru¬ 
cially  important  for  the  enterprise  to  continue  [11,  8,  12,  13].  Because  of  the 
difficulties  of  quantum  coherence,  the  first  quantum  computer  comprised  only 
two  qubits. 

An  historical  starting  point  that  led  to  quantum  computing  was  reversible 
computing  [14].  Since  microscopic  physics  is  reversible4,  it  is  believed  that 
quantum  mechanical  algorithms  must  be  too.5  Reversible  algorithms  for  sim¬ 
ulating  physics  on  a  quantum  device  can  serve  as  a  guide  for  constructing  the 
device.  The  common  assumption  is  the  quantum  mechanical  device  itself  un¬ 
dergoes  unitary  (and  therefore  reversible)  evolution  as  it  transitions  through  its 
“computation”.6 

For  any  reversible  computation,  one  can  describe  the  algorithm  by  specify¬ 
ing  a  unitary  evolution  operator,  formally  written  as  e'^r/h,  acting  on  the  sys¬ 
tem  wavefunction,  1$),  which  constitutes  the  state  of  the  quantum  computer’s 
“memory”.  With  N  qubits,  the  quantum  state  |4>)  resides  in  an  exponentially 
large  Hilbert  space  with  2N  dimensions.  A  new  quantum  state,  | '!>'),  is  generated 
by  application  of  a  unitary  matrix  of  size  2N  x  2N  as  follows 


|$')  (1) 

By  repeated  application  of  e'Hrth  an  ordered  sequence  of  states  is  generated 
and  each  one  is  given  a  unique  time  label.  If  the  first  state  is  labeled  by  t  then 
the  next  one  is  labeled  by  t  +  r,  and  the  next  by  t  +  2 r,  and  so  forth.  With  this 

2 Cory  et  al.  have  employed  the  quantum  number  mz  of  a  nuclear  spin  of  an  atom  in 
a  molecule  of  a  liquid  placed  in  a  strong  external  magnetic  field  to  encode  a  single  qubit 
and  they  used  nuclear  magnetic  resonance  to  control  its  state  and  interaction  with  qubits  in 
neighboring  atoms  within  the  same  molecule  [10]. 

3In  lattice-gas  quantum  computation,  quantum  parallism  is  used  to  allow  for  simultaneous 
multiple  collision  possibilities  at  each  site  of  the  lattice.  This  allows  for  a  reduction  in  the 
viscosity  of  the  fluid  that  improves  the  computational  efficiency.  Yet  the  computational  effi¬ 
ciency  is  also  due  to  the  continuous  phase  of  the  qubit,  |<?)  =  cos0|l)  +  sin0|O),  which  we  use 
to  represent  the  probability  of  finding  a  particle,  fa  =  cos2  0. 

4  Provided  photons  do  not  escape  to  infinity. 

5 In  this  paper  we  use  irreversibility,  for  practical  purposes,  in  part  of  the  quantum  me¬ 
chanical  algorithm. 

6By  restricting  oneself  to  reversible  algorithms,  in  principle  heat  production  may  be  avoided 
altogether  [15,  16]. 
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understanding  we  write  (1)  as 

|*(t  +  r))=eiAr/ft|tf(t)).  (2) 

In  this  way  the  computational  time  advances  incrementally  in  unit  steps  of  du¬ 
ration  t.  Of  course  the  state  of  the  quantum  computer  exists  at  all  intermediate 
times,  say  at  t+%,  but  for  our  purposes  we  need  only  use  the  state  at  intervals  of 
the  time  step  r.  The  quantum  computer’s  evolution  is  invertible  by  application 
of  the  adjoint  of  the  evolution  operator 

|  $(t  -  t))  =  e~iftr/h  |  $(*)).  (3) 

This  computational  picture  is  consistent  with  the  Heisenberg  picture  of  quantum 
mechanics.  For  any  reversible  algorithm  chosen,  the  task  is  to  map  the  algorithm 
onto  the  dynamical  evolution  of  interacting  qubits  of  the  physical  device. 

3  Quantum  Lattice  Gas 

Lattice-gas  quantum  computation  uses  the  superposition  of  multiple  qubit  states 
within  a  small  spatial  region  of  size  £  only  for  a  short  amount  of  time  on  the 
order  of  the  duration  of  a  single  time  step,  r.  A  lattice-gas  quantum  computer 
has  qubits  arranged  in  a  lattice-based  array,  with  a  small  group  of  qubits  at  each 
site.  Each  site  of  the  lattice  can  be  thought  of  as  a  small  quantum  computer  and 
all  the  quantum  computers  are  connected  in  a  lattice  network.  The  quantum 
lattice  gas’  evolution  can  be  formally  expressed  as  a  special  case  of  (2)  as  follows 

|  ^(fi,...,xv;t  +  r))  =  SC  |  ^(xi,...,xv;*)>,  (4) 

where  S  is  a  unitary  streaming  matrix  and  C  is  a  unitary  collision  matrix  and 
where  we  have  explicitly  labeled  the  wavefunction’s  dependence  on  all  the  co¬ 
ordinates  of  the  lattice.  The  operator  C  causes  mixing  of  the  outgoing  collision 
configuration  at  each  site  of  the  lattice,  locally  entangling  the  qubit  states  within 
a  lattice  cell  of  radius  size,  i.  The  operator  S  causes  qubits  to  move  from  one 
site  to  the  next,  by  exchanges  between  nearest  neighboring  sites  (it  is  identical 
to  its  classical  counterpart).  Each  qubit  moves  with  unit  speed,  c  =  £,  along 
one  of  the  lattice  directions,  ea.  Hence,  in  a  completely  coherent  quantum  com¬ 
putation,  the  application  of  S  causes  global  entanglement  of  the  all  the  qubit 
states7.  It  remains  an  intractable  problem  to  theoretically  analyze  the  dynamics 
of  a  quantum  computer  with  many  qubits  because  of  the  exponentially  large  size 
of  the  Hilbert  space  in  which  the  entanglement  occurs.  And  to  make  matters 
worse,  even  if  a  quantum  computer  was  constructed  with  a  large  number  of 
qubits,  its  wavefunction  would  decohere  by  uncontrolled  entanglement  with  the 
external  world  and  we  know  of  no  way  to  mitigate  against  this.  So  constructing 

7Mathematically  speaking,  this  is  because  both  5  and  C  cannot  be  simultaneously  diago¬ 
nalized 
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a  large  coherent  quantum  computer  is  difficult,  if  not  all  together  impossible, 
and  predicting  its  behavior  by  analytical  means  is  intractable. 

So  what  can  be  done  about  this?  What  I  would  like  to  consider  is  a  simplifi¬ 
cation  that  will  sidestep  these  obstacles  and  give  us  two  important  advantages: 
(1)  a  simple  way  to  use  a  quantum  computer  with  a  large  number  of  qubits; 
and  (2)  a  way  to  analysis  of  its  behavior.  In  lattice-gas  quantum  computation 
complete  coherence  of  the  wavefunction  is  not  needed  for  the  algorithm  to  work. 
In  fact,  we  assume  entanglement  of  qubit  states  is  only  among  small  clusters  of 
qubits  in  a  localized  nearby  neighborhood,  so  independent  quantum  operations 
are  done  in  a  classically  parallel  fashion  on  all  sites  simultaneously.  This  is  the 
collision  step. 

In  a  deterministic  classical  lattice  gas,  the  collision  operator  is  a  permutation 
matrix  with  components  being  either  zero  or  one.  In  a  probabilistic  classical 
lattice  gas,  the  collision  operator  is  a  transition  matrix  with  real  valued  com¬ 
ponents.  In  contrast,  in  a  quantum  lattice  gas,  the  collision  operator  can  be  a 
unitary  matrix  with  complex  components.  The  collision  process  is  in  general 
irreversible  because  a  projection  of  the  quantum  computer’s  wavefunction  into 
a  tensor  product  state  over  the  qubits  is  periodically  made  causing  the  wave- 
function  to  partially  collapse.  Hence,  application  of  S  does  not  cause  any  global 
entanglement. 

The  quantum  lattice  gas  presented  here  should  not  be  confused  with  previ¬ 
ous  quantum  lattice  gas  models  by  Succi  [17],  Boghosian  [18],  or  Yepez  [19]  for 
simulating  quantum  mechanical  systems.  Despite  some  similarities,  the  type  of 
quantum  lattice  gas  treated  in  this  paper  is  a  direct  generalization  of  a  classi¬ 
cal  lattice  gas  with  quantum  bits  replacing  classical  bits.  In  fact,  if  orthogonal 
permutation  matrices  with  0  and  ±1  components  are  used  for  the  collision  pro¬ 
cess,  in  the  limit  of  complete  collapse  of  the  lattice-gas  quantum  computer’s 
wavefunction,  the  quantum  lattice  gas  exactly  reduces  to  a  classical  lattice  gas. 
This  particular  feature  distinguishes  the  quantum  lattice  gas  for  fluid  simulation 
from  the  quantum  lattice  gases  for  quantum  mechanical  simulation. 

4  Preliminaries 

Consider  a  lattice-gas  quantum  computer  with  the  following  properties: 

•  V  is  the  number  of  lattice  sites 

•  B  is  the  number  of  qubits  per  site  (and  the  number  of  nearest  neighbors) 

•  N  =  VB  is  the  total  number  of  qubits 

•  2N  is  the  size  of  the  full  Hilbert  space 

•  2b  is  the  size  of  the  on-site  submanifold,  denoted  H 

•  B  is  the  size  of  the  reduced  on-site  submanifold,  denoted  B 
We  will  use  the  following  convention  for  indices: 
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Constants 

Names 

e 

length  unit 

r 

time  unit 

m 

mass  unit 

c 

velocity  unit  (f) 

D 

spatial  dimension 

B 

lattice  coordination  number 

€-ai 

unit  lattice  vectors 

a 

directional  index  (1,2,. . . ,  B) 

i,j,k,l 

spatial  indices 

Table  1:  Model  Constants 


Table  2:  Wavefunction  Symbols 


Symbol 

Size  of  Manifold 

Description 

2^ 

Total  system  wavefunction 

2b 

On-site  ket 

U) 

B 

Partially  collapsed  on-site  amplitudes 

Q 

2 

Qubit  ket 

•  Small  roman  letters  (a,  b,  c)  for  the  5-space  dimensions,  a  6  {1, . . . ,  B} 

•  Greek  letters  (a,  (3,  7)  for  the  %-space  dimensions,  a  6  {0, . . . ,  2B  —  1} 

•  Middle  roman  letters  ( i ,  j,  k)  for  the  spatial  dimensions,  i  €  {1, . . . ,  D} 

The  full  Hilbert  space  of  size  2BV  is  partitioned  into  V  independent  quantum 
manifolds  of  dimension  2s,  as  depicted  in  Fig.l.  Quantum  superposition  of 
states  occurs  only  within  each  2s-dimensional  subspace,  denoted  H.  A  general 
on-site  ket  defined  over  the  basis  states  of  Ti  is  the  following 


o=0 


(t\ 


(5) 


The  ket  |  ip)  is  specified  by  2s  complex  amplitudes,  denoted  ipo, . . . ,  V>2B- 1-  The 
quantum  computer’s  total  wavefunction  is  formed  as  a  tensor  product  over  all 
the  'H-manifolds 

v 

|  V(xi,...,xv;t))  =  (^)  |  rp(x,t)).  (6) 

X=1 

The  collision  operator,  C,  is  blocked  over  all  the  "H-manifolds.  That  is,  the 
collision  matrix  is  block  diagonal  with  V  blocks  each  of  size  2B  x  2s,  and 
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Lattice  Network 


*  J* 


Single  Lattice  Site 


Figure  1:  An  array  of  small  quantum  computers  (the  quantum  computers  are  depicted  as 
circles)  arranged  in  a  2-dimensional  triangular  lattice  (B=6).  The  large  circle  on  the  right  is 
an  expanded  view  of  a  single  quantum  computer,  which  is  one  site  of  the  lattice.  It  depicts 
the  on-site  submanifold,  'H.  Each  quantum  computer  at  a  lattice  node  has  6  qubits  so  the 
on-site  ket  |  V>)  resides  in  a  64-dimensional  Hilbert  space.  Each  node  is  coupled  to  its  6  nearest 
neighboring  quantum  computers  by  a  mechanism  allowing  for  the  exchange  of  a  single  qubit. 


therefore  can  be  written  as  a  tensor  product 

v 

C  =  <g)U.  (7) 

x—\ 

The  on-site  collision  matrix,  U,  is  unitary  and  acts  on  the  on-site  ket 

\ip'(x,t))  =  U\rp(x,t)).  (8) 

The  prime  on  the  L.H.S.  of  (8)  indicates  that  the  ket  is  an  outgoing  collisional 
state. 


5  Unitary  Collision  Matrix 


Let  Qa  be  matrices  representing  the  conserved  quantities  in  the  single  speed 
quantum  lattice  gas,  Qa  =  ( Qo,Qi )  where  i  is  an  index  over  the  independent 
spatial  coordinates.  A  fundamental  property  of  a  quantum  lattice  gas  is  that 
the  mass  density  and  the  momentum  density  can  be  written  as  follows 


P  =  (ip  I  Qo  | 
pvi  =  {ip  |  Qi  |  ik), 


that  is,  where  the  component  of  Qa  are 


B 

(Qo)pv  rnS^i,  ^  ^ 

a=l 


(9) 

(10) 


(11) 
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and  where  the  component  of  Qi  are 

B 

{Qi'itiv  =  TflcS^  ^  '  bflnCal ■  (12) 

a=l 


The  b^a  denotes  the  ath-bit  of  the  fith  ket.  The  ea  here  denote  the  unit  lattice 
vectors  where  a  =  1  The  components  of  the  2B  x  2s  qubit  number 

operator  are  defined  by 


In  terms  of  (13),  the  operators  for  mass  and  momentum  are 

(13) 

■o> 

o 

II 

(14) 

and 

B 

Qi  —  TTIC  ^  ^  eaj7la . 

a=l 

(15) 

In  terms  of  (13)  the  invariant  quantities  are  simply  expressed  as 
matrix  elements 

the  following 

B 

P  =  m &  1  i 

a=  1 

(16) 

B 

pVi  =  ^2mceai(ip  I  nQ  |  ip). 

a=l 

(17) 

The  matrix  elements  (9)  and  (10)  must  remain  constant  after  each  time  step 
iteration 

(ip{t  +  r)  |  Qa  |  i>{t  +  r))  =  (ip(t)  |  Qa  |  ip{t)).  (18) 

Since  |  xp(t  +  r))  =  U  |  this  implies  that 

V'QaU  =  Qa, 

(19) 

which  is  just  the  commutator 

[£,gtt]  =  o. 

(20) 

The  matrices  Qa  must  commute  with  U. 

Let  g  denote  the  generator  of  U 

U  =  ei£\ 

(21) 

where  e  is  an  “Euler  angle”.  Consider  a  “rotation”  through  an 
angle  e  so  that  U  can  be  Taylor  expanded  to  first  order  as 

infinitesimal 

U  =  1  +  ieg. 

(22) 
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The  unitary  condition,  U W  =  1,  implies  that  the  generator  is  hermitian 

g-g'  =  0  +  O(e2).  (23) 

Prom  (19),  we  see  that  mass  and  momentum  conservation  is  ensured  provided 

Qag-g'Qa  =  0  +  O(e2).  (24) 

The  solution  of  the  set  of  linear  equations  (23)  and  (24)  give  the  Lie  algebra  for 
the  unitary  group.  Therefore,  the  mass  density  (9)  and  the  momentum  density 
(10)  are  conserved  when  each  equivalence  class  block  of  the  collision  operator  is 
an  element  of  the  unitary  group  U ( n )  where  n  is  the  size  of  the  equivalence  class 
of  the  incoming  local  configuration.  This  is  an  important  feature  of  a  quantum 
lattice  gas.  Since  any  member  of  the  unitary  group  can  be  used,  the  quantum 
lattice  gas  is  algorithmically  robust. 


Figure  2:  The  equivalence  classes  for  the  quantum  FHP  lattice  gas. 


An  equivalence  class  is  defined  as  a  set  of  basis  states  that  correspond  to  par¬ 
ticle  configurations  with  the  same  mass  and  momentum.  The  unitary  collision 
operator,  U,  acting  on  the  2s  dimensional  7f-manifold  itself  is  block  diagonal 
over  the  equivalence  classes.  For  example,  there  are  two  equivalence  classes  for 
the  FHP  lattice  gas  [2],  see  Fig. 2.  The  first  equivalence  class  is  comprised  of  the 
following  two-body  kets 

1 9)  =  1001001) 

1 18)  =  1010010) 

|36)  =  j  1001 00) 
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with  mass,  m  =  2,  and  zero  momentum,  p  =  0.  A  general  ket  in  this  mass- 
momentum  sector  of  the  on-site  manifold  is  a  linear  combination  of  these 

a\  100100)  +  /?|010010)  +  7)001001),  (25) 

where  a,  /?,  and  7  are  complex  numbers.  The  second  equivalence  class  is  com¬ 
prised  of  the  following  three-body  kets 

|21)  =  1010101) 

|42)  =  1101010) 

with  mass,  m  —  3,  and  zero  momentum,  p  =  0.  A  general  ket  in  this  mass- 
momentum  sector  is  a  linear  combination  of  these 

At|101010)  +  ^JOIOIOI).  (26) 

So  U  for  a  two-dimensional  quantum  lattice  gas  on  a  triangular  lattice  has  two 
blocks,  a  17(3)  block  for  mixing  the  2-body  configurations  and  a  U{ 2)  block  for 
mixing  the  3-body  configurations. 

For  the  triangular  quantum  lattice  gas,  we  have 

(**')- «"  (  e“ts,n’'  )  (M  ,  (27) 

\V>42  /  \-e  4  sin r?  e  Hcosj?y  \1p42  J 

where  zero  momentum  three-body  configurations  are  mixed  by  a  unitary  matrix, 
U( 2)  =  1/(1)  ®  SU (2),  which  in  general  has  four  free  parameters.  The  zero 
momentum  two-body  configurations  are  mixed  by  a  unitary  matrix,  U( 3)  = 
U(  1)  ®  517(3),  which  in  general  has  nine  free  paramters8 

Uis  =  eieSU(3)  U18  .  (28) 

\^36  /  V^36/ 


6  Partial  Collapse  of  Post- Collision  Ket  |  ip') 


To  avoid  causing  any  global  entanglement  as  induced  by  streaming,  we  project 
the  post-collision  ket  |  ip1}  which  resides  in  the  2s-dimensional  'H  manifold 
onto  a  smaller  B-dimensional  submanifold,  B.  This  is  done  using  a  projection 
operator,  denoted  T,  as  follows 


I  “/>  =  f  |  it) 


(29) 


The  operator  f  causes  a  partial  collapse  of  the  locally  entangled  on-site  state  | ip) 
resulting  in  a  nonentangled  state  |u;)  residing  in  a  smaller  manifold  (a  mapping 
8We  do  not  write  out  the  SU(3)  matrix  in  component  form  because  it  is  too  complicated. 
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Dimensions 


B"  Dimensions 


H 

xj)3 


000000> 

000001> 

000010> 

000011> 
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111110> 


Projection 


Injection 


Figure  3:  Two-dimensional  quantum  lattice  gas  on  a  triangular  lattice.  The  lattice  coordi¬ 
nation  number  is  B  =  6.  There  are  26  =  64  amplitudes  in  the  'H-manifold  and  6  amplitudes  in 
the  0-manifold.  Injection  maps  the  6  on-site  amplitudes  in  the  0-manifold  into  the  larger 
W-manifold.  The  inverse  process,  projection,  maps  the  64  amplitudes  ipa  in  the  'H-manifold 
onto  6  amplitudes  ua  in  the  0-manifold.  The  projection  is  a  measurement  process  that  causes 
a  partial  collapse  of  the  on-site  wavefunction  |  ip(x,  t)).  The  partially  collapsed  wavefunction 
is  |  u>(x,  t)). 


from  64  dimensions  down  to  6,  see  Fig.3).  Thus  |w)  in  (29)  may  be  termed  the 
collapsed  post-collisional  on-site  ket.  By  construction,  the  action  of  F  fixes  the 
phase,  6a,  of  the  on-site  qubits  |  qa)  according  to  the  following  recipe 

|  qa)  =  cos 0a  |  1)  +  sin0o  |  0),  (30) 

where  u>a  =  cos  0a .  That  is, 

I  9o)  =  wa  I  1)  +  \/l  -  I  0).  (31) 

After  the  collapse  of  the  ket  |  ip)  the  single-particle  occupation  probability,  de¬ 
noted  fa,  is  a  well-defined  quantity.  It  is  the  probability  of  finding  a  particle  at 
coordinate  (x,  t)  with  momentum  mcea 

fa{x,  t)  =  u*a(x,  t)u>a(x,  t).  (32) 

Using  the  single  qubit  number  operator  h  =  ^  ^  ^  ^ ,  (32)  can  be  written  in 
terms  of  the  qubit  ket  as 

fa(x,t)  =  ( qa{x,t )  |  h  |  qa(x,t)).  (33) 

Furthermore,  in  (16)  and  (17),  the  matrix  element  (ip  \  ha  \  ip)  also  gives  the 
probability  of  particle  occupancy.  So  for  a  quantum  lattice  gas,  fa  can  also  be 
expressed  as  the  matrix  element  of  the  multiple  qubit  number  operator9 

_ fa  =  I  na  I  ip).  (34) 

9In  a  classical  lattice  gas  the  single-particle  occupation  probability  is  obtained  by  ensemble 
averaging  over  the  number  variables,  /„  =  (na),  where  n„  =  0  or  1. 
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Inserting  the  expression  for  the  outgoing  collisional  state  (8)  into  the  R.H.S. 
of  (29),  we  have 

|  u>')  =  TU  |  tp).  (35) 

In  §7  we  will  use  a  nonlinear  function  for  the  projection  operator.  Our  goal 
is  to  retain  as  much  quantum  information  as  possible  while  allowing  (35)  to 
reduce  to  the  collision  equation  of  classical  lattice  gas  transport  when  U  is  a 
real-valued  permutation  matrix.  This  is  accomplished  by  projecting  down  from 
the  "H-manifold  containing  2s  complex  amplitudes  to  the  6-manifold  with  only 
B  complex  amplitudes.  Each  uia  (or  associated  qubit  |  qa))  is  “attached”  to 
one  of  the  lattice  directions.  The  reason  for  this  reduction  of  the  quantum 
information  is  the  following.  By  using  only  B  complex  amplitudes,  one  for 
each  direction,  we  can  straightforwardly  write  down  a  quasi-classical  streaming 
equation  in  analogy  to  the  streaming  equation  of  a  classical  lattice  gas 

|  uj(x  +  eCa,t  +  r))  =|  u>'(x,t)),  (36) 

where  Cabi  =  eai?iah.  Inserting  (35)  into  (36),  we  have 

\uj(x  +  e£a,t  +  r))  =TU\1>(x,t)).  (37) 

The  only  task  remaining  to  complete  the  analogy  to  classical  lattice  gas  dy¬ 
namics  is  to  rewrite  the  R.H.S.  of  (37)  solely  in  terms  of  the  w’s.  This  can  be 
done  by  injecting  the  on-site  collapsed  ket  |  w)  residing  in  the  sub-manifold  B 
up  into  the  larger  on-site  manifold  V.  (see  Fig.3).  This  process  can  be  expressed 
by  application  of  an  injection  operator,  I,  as  follows 

|  Tp(x,t))=I  \u(x,t)).  (38) 

A  straightforward  way  to  accomplish  the  injection  is  to  take  the  tensor  product 
over  the  on-site  qubits 


B 

I  i>(x,t))  =  (g)  |  qa(x,t)).  (39) 

a— 1 

This  is  a  nonlinear  operation.10  Let  us  revisit  the  example  a  two-dimensional 
quantum  lattice  gas  on  a  triangular  lattice  ( B  =  6),  a  generalization  of  the 
classical  FHP  lattice  gas  [2],  Fig.3  illustrates  projection  from  the  64-dimensional 
Tf-manifold  down  to  the  6-dimensional  B-manifold  and  illustrates  injection  from 
the  6-manifold  up  to  the  -manifold.  The  6  on-site  amplitudes  u>0  (or  the 
associated  6  on-site  qubits  |  qa))  generated  by  the  projection  can  be  streamed 
in  a  classical  fashion.  After  streaming  to  their  new  sites,  each  quantum  computer 
has  a  new  incoming  configuration  of  the  u>a  amplitudes.  Before  this  configuration 
can  be  collided,  they  must  be  injected  up  to  the  larger  64-dimensional  manifold 
where  the  collision  process  is  well  defined. 

10 A  non-square  linear  matrix  could  also  be  used,  but  it  is  difficult  to  find  an  appropriate 
matrix  even  though  it  can  be  shown  that  one  exists. 
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(40) 


Inserting  (38)  into  (35)  gives 

\u,'(x,t))  =  rui\oJ(x,t)). 

Using  the  fact  that  the  projection  of  the  injection  is  the  identity  operation: 
|  ui)  =  VI  |  w),  we  write  (40)  in  a  form  analogous  to  the  classical  lattice  gas 
collision  equation 

I  u/(f,  t))  =|  u(x,  t))  +  [f Ui  I  w(£,  t))  -  ti  I  w(x,  *)>]  .  (41) 

Finally,  we  arrive  at  the  quantum  lattice  gas  microscopic  transport  equation  by 
inserting  (41)  into  (36) 

|  ui{x  +  £Ca,t  +  t))  =|  u>(x,i))+  |  Sl{x,t)),  (42) 

where  the  quantum  lattice  gas  collision  operator  is  defined  as 

|  fi(x,  t))  =  ft//  |  u{x,  t))  -  ti  |  w(x,  t)).  (43) 

In  component  form,  (42)  is 

u>a(x  +  tea,t  +  t)  =  uia(x,t)  +  Qa  (w*(x,t))  .  (44) 

Equation  (44)  is  identical  in  form  to  the  classical  lattice  gas  transport  equation 
where  the  occupation  variable,  na  =  0  or  1,  is  replaced  by  a  complex  amplitude, 
0  <  |wa|  <  1,  that  continuously  encodes  the  square  root  of  the  probability  for 
particle  occupancy.  Hence  (44)  is  a  much  more  useful  expression  of  the  quantum 
lattice  gas  dynamics  than  (4)  is. 

7  The  Projection  Operator 

We  can  write  an  analytical  expression  for  the  projection  operator  where  the 
amplitudes  oja  a  nonlinear  function  of  the  amplitudes  ipa 


Note  that  baa  =  0  or  1  is  the  Boolean  value  of  the  a*’’  bit  of  the  ath  ket  in  the 
number  representation.  Let  m  and  p*  be  the  operators  for  mass  and  momentum 
in  the  S-space.  Then  the  matrix  element  for  the  mass  density  is 

p=(w|m|w),  (46) 

where  mab  —  rn6at,,  and  the  matrix  element  for  the  momentum  density  is 

pvi  =  {w\pi\  w),  (47) 
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where  ( Pi)ab  =  mceai8ab ■  Qa  and  Qi  were  defined  in  §5  to  be  the  mass  and 
momentum  operators  in  V'-space.  Here  we  have  mass  and  momentum  operators 
in  w-space.  The  matrix  element  (9)  defines  the  mass  density  as  p  =  {ip  \  Q0  | 
ip),  where  (Q0)ap  =  mY^=\baa5ap,  and  the  matrix  element  (10)  defines  the 
momentum  density  as  pvt  =  (ip  |  Q;  |  ip),  where  {Qi)ap  =  mc^2^=1baaeaiSab. 
Equating  (46)  with  (9)  and  equating  (47)  with  (10)  gives  us  a  way  to  check  the 
projection  operator  (45). 

This  is  done  as  follows 


B 


(w  |  m  \  u>)  =  m  ^  |  u>a  |2 

a— 1 

(48) 

B  2b 

=  ^2  I  ipa  I2  baa 

0=1 

(49) 

2b  (  B  \ 

=  E  i  ^  i2  mE6“° )  ’ 

a=l  \  a= 1  / 

(50) 

where  we  used  the  square  of  (45)  on  the  second  line  of  the  derivation.  Therefore, 
we  have 

{u)\m\u)  =  (ip\Qo\  ip),  (51) 

where  the  mass  operator  in  i/'-space 

B 

{Qo)aP  —  ^  '  baa5ap 

(52) 

a=l 


is  identical  to  Q0  defined  in  (9).  So  the  projection  operator  (45)  conserves  mass. 
We  continue  the  consistency  check  by  rewriting  (47) 


(w|pj|w)  =  mc]T|wa|2ea; 

a=l 
B  2b 

=  me  EE  I  "*Pa  |  bQaea 


a=  l  a=l 


—  ^  ]  |  Ipa  |  |  TTIC  ^  ^  baaeai  |  , 


a=l 


(53) 

(54) 

(55) 


where  again  we  used  the  square  of  (45)  on  the  second  line  of  the  derivation. 
Therefore,  we  have 

(w  |  ^  |  w)  =  {ip  |  Qi  |  ip),  (56) 

where  the  momentum  operator  in  ip- space 


{Qi)a0  —  TTIC  'y  KfbaaCajSap, 


(57) 
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is  identical  to  Qi  defined  in  (10).  So  the  projection  operator  (45)  conserves 
momentum  as  well  as  mass. 


8  Equilibrium  Ansatz 

The  u!a  amplitudes  in  B-space  can  be  ordered  in  powers  of  e  as  follows 

|  oj)  =|  u^°))  +  e  |  u/1))  +  0(e2),  (58) 

where  |  k/0))  denotes  the  equilibrium  ket  and  where  the  ket  |  u>M)  is  the  first 
order  correction  from  equilibrium.  The  condition  equilibrium  is  that  |  a/0)) 
satisfies  the  following  identity 

|  u/°))  =  ff// |  w(0)).  (59) 

Note  that  the  associated  equilibrium  ket  in  7i- space,  |  tp^),  follows  from  (58) 
by  injection  |  ip)  =  7  |  ui) 

|V0=l  ^0))  +  e  \  iPw)  +  0(e2).  (60) 

So  we  can  also  write  (59)  as  follows 

|  V>(0))  =  U  |  V-(0)).  (61) 

It  is  clear  that  |  tp^)  is  an  eigenvector  of  U  with  unity  eigenvalue.  Using  (59), 
we  immediately  see  that  the  collision  operator  (43)  vanishes  at  equilibrium 

|  fi(w<0)))  =  Wl  |  w<°>)  -  tl  |  w(0))  =  0.  (62) 

Equations  (58)  and  (59)  constitute  the  essential  ansatz  that  will  allow  us  to 
perform  a  Chapman-Enskog  analysis  of  the  quantum  lattice  gas.  It  is  possible  to 
analytically  solve  (59)  for  |  ui^).  Knowing  the  form  of  |  w<°>)  ,  we  can  predict  the 
hydrodynamics  equations  of  the  quantum  lattice  gas  at  the  macroscopic  scale. 
In  the  Chapman-Enskog  analysis,  we  expand  the  collision  operator,  |  0) ,  about 
this  equilibrium  ket  |  u/°)).  In  so  doing,  the  Jacobian  of  the  collision  operator 
is  computed  as  a  first  order  correction  and  is  evaluated  at  |  oj)  =|  a/0)).  The 
transport  coefficients  for  the  mass  diffusion,  shear  viscosity,  and  bulk  viscosity 
depend  on  the  value  of  this  first  order  correction,  and  this  in  turn  depends  on 
the  value  of  |  w^0').  Hence,  one  must  determine  the  equilibrium  amplitudes  in 
order  to  compute  the  value  of  the  transport  coefficients. 

The  a  single  particle  occupancy  probability  fa  =  |  h  \ 

q^)  =  {%p  |  ha  |  ip)  has  the  functional  form 

~  gap+0ea-p+~fE  _|_  2  ’ 

where  the  argument  of  the  exponential  is  a  linear  combination  of  the  conserved 
scalar  quantities:  (1)  the  mass  p;  (2)  the  momentum  component  ea  ■  p  along 
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Figure  4:  Theory  versus  simulation  comparison  of  the  velocity  dependence  of  the  single¬ 
particle  distribution  function  in  the  non-Galilean  parameterization:  fa  =  u>*o>a  =  d  +  dDea  ■ 
v  +  gdD(D/2  +  l)Qa  :  vv.  FHP  simulation  data  is  overplotted  on  this  predicted  mesoscopic 
distribution  function.  Plots  (a)  and  (b)  are  for  background  densities  of  d  =  .20  and  d  =  0.25, 
respectively.  A  velocity  shift  is  imparted  along  the  x-axis;  that  is,  along  the  fi  direction 
indicated  in  the  figure.  Data  was  collected  from  a  128  X  128  classical  FHP  simulation  (crosses) 
and  was  coarse-grained  averaged  over  1600  time  steps  from  time  step  t  —  400  to  t  =  2000. 
Data  was  also  collected  from  a  smaller  32  X  32  quantum  FHP  simulation  (circles)  and  were 
measured  at  a  single  time  step  at  t  =  200. 


the  lattice  direction  ea;  and  (3)  the  energy  E  at  a  lattice  site.  The  real  valued 
coefficients  a,  /3,  and  7  are  free  parameters  that  are  fixed  by  the  non-Galilean 
parameterization  given  in  Appendix  A.  The  reason  for  the  form  of  (63)  is  the 
collision  matrix  U  is  unitary  and  so  the  collisions  obey  detailed  balance. 

In  the  quantum  limit,  where  the  quantum  lattice  gas  becomes  a  fully  coherent 
quantum  system  that  undergoes  unitary  evolution,  we  expect  the  equilibrium 
probability  for  the  particle  occupancy  to  have  the  form  of  (63).  In  the  oppo¬ 
site  limit,  the  classical  limit,  where  there  is  a  complete  collapse  of  the  ket  |  'P) 
everywhere,  the  quantum  lattice  gas  reduces  to  a  classical  lattice  gas  system. 
In  the  classical  limit  too,  the  particle  occupancies  are  described  by  (63).  Our 
quantum  lattice  gas  dynamics  is  somewhere  midway  between  a  fully  quantum 
system  and  a  fully  classical  system.  Since  both  ends  of  the  spectrum  are  de¬ 
scribed  by  (63),  it  is  not  altogether  unexpected  that  the  middle  regime  is  too. 
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However,  this  is  not  obvious.  At  first  glance,  it  appears  that  the  destructive  act 
of  projecting  from  the  %-manifold  onto  the  B-manifold  might  destroy  the  form 
of  the  equilibrium  distribution.  But  it  does  not  (see  Fig.4);  the  distribution 
(63)  describes  the  particle  occupancy  of  our  quantum  lattice  gas  [3].  Taking 

u>a  =  \J  {q^  |  fi  |  q'a'1),  the  result  (see  Appendix  A)  is  the  following11 

I  «*>)  =  Vd  1)  +  ^ v ■  |  e)  +  9D^c2  -vv  :|  Q)j  ,  (64) 


where 

/!\ 
id-  ; 

i 

Vi/ 

ID-^ID, 

and 

|Q)  =  (££-i)|l>. 

Taylor  expanding  the  collision  operator  gives 


(65) 


(66) 

(67) 


fia  (u>)  =  Ha  +  £ 


dQ.a  (w) 


dojb 


J=u>m 


4l)  +  o(e2). 


(68) 


Using  the  equilibrium  condition  (62),  the  first  term  on  the  R.H.S.  vanishes  and 
we  are  left  with 

na(oJ)  =  £jabJb1)  +  0(e2),  (69) 

where  the  Jacobian  of  the  projection  operator  is  defined  as 


anQ(w) 

dab  — 


dujb 


u»=u>(°) 


Equation  (69)  can  be  written  in  vector  form  as 

|  Cl  (w))  =  ej  |  u>w)  +  0{e2). 


(70) 


(71) 


9  Boltzmann  Equation  for  the  Quantum  Lattice 
Gas 

The  quantum  lattice  gas  microscopic  transport  equation  (42)  is 

|  u>(xl  +  e£C,  t  +  £2t))  =|  w(x,  t))+  |  f2(x,  t)). 

11In  (64),  we  have  used  the  approximation  vd  +  x  —  1  +  ^ ,  which  holds  for  small  x. 
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Near  equilbrium,  each  application  of  the  projection  operator  causes  only  small 
changes  in  the  phase  of  |  uj)  allowing  us  to  Taylor  expand  the  L.H.S.  of  (42). 
Taylor  expanding  (42)  to  second  order  in  e  gives  a  Boltzmann  equation 

(a  c-2/2  A  A  \  1 

ecCidi  +  -^—CiCjdidj  j  |  uj(x,t))  +  0(e3)  =  -  \  Q(x,t)). 

J  (72) 

In  analogy  with  a  classical  lattice,  we  impose  two  constraints  on  the  collision 
operator  regarding  the  isometries  of  the  lattice.  The  first  constraint  is 

|  m  |  ft)  =  0,  (73) 


and  this  will  be  needed  when  we  take  the  zeroth  moment  of  (72).  This  constraint 
enforces  mass  conservation.  The  second  constraint  to  enforce  momentum  con¬ 
servation  is 

<«|  ft  |  ft)  =  0,  (74) 

and  this  will  be  needed  when  we  take  the  first  moment  of  (72).  Inserting  the 
e-expansion  of  |  ua),  (58),  and  the  e-expansion  of  |  fiQ),  (71),  into  the  quantum 
lattice  gas  transport  equation,  (72),  and  keeping  terms  up  to  second  order  in  e 
gives 


e2at  I  I  u »W)+— Aft  I  <* >w)+^-cicjaidj  \  «<°>)  =  -j 

T  T  IT  T 

Now  equating  the  order-e  terms,  we  get 


w(1  ))+P(e3). 
(75) 


J  |  w^))  =  Cidi  |  w(0)). 


(76) 


Then  inverting  the  kinetic  part  of  the  Jacobian  matrix  gives 

|  w*1*}  =  ej~l  CA  |  w<°)).  (77) 

So  the  first  order  correction  ket  is  equated  to  the  gradient  of  the  equilibrium 
ket.  Inserting  this  result  back  into  (58)  we  have 

|  u,)  =|  «<°>)  +  ej-'tidi  |  w{°>)  +  0(e2),  (78) 

Inserting  (64)  into  (78),  we  have  the  subsonic  e-expansion  of  |  u>)  good  to  first 
order  in  e 

\u)  =  >/d(\l)+  °v-  |  e)  +  9D{V  +  2).vv  :|  Q)  +  ^  J~lCC  :  Vu  |  1))  +0(e2). 

(79) 


10  Hydrodynamic  Equations 

ZEROTH  MOMENT  EQUATION: 
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Now  we  can  write  the  zeroth  moment  of  the  Boltzmann  equation  (72)  by  left 
multiplying  by  (o>|m  as  follows 


.  £2^2  A  „  i 

e2{ui  |  m  |  dtoS)  +e(w  |  mcCi  |  d;u>)  +  — — (ui  |  mtitj  \  didju )  =  —  (u>  |  rh  |  fl). 

2t  t 

(80) 

The  adjoint  of  (80)  is  the  following 


a  £-2^2  a  A  2 

e2(dtu>  |  rh  |  w)  +  e{diu>  |  Lithe  |  a;)  +  — — ( didjui  |  tiCjih  |  u)  =  -  (ft  |  m  |  u>), 

2r  t 

(81) 

and  the  R.H.S.  will  also  vanishes  because  of  constraint  (73).  Because  m  and  £ 
commute,  adding  (80)  and  (81)  gives  us  the  following  zeroth  moment  equation 

A  £-2^2  /  A  A  A  A  \ 

e2dt(oj  |  m  |  u>)+£di(w  |  mcCi  |  w)+— —  |  mtitj  |  d{ djU>)  +  {d%djU>  |  titjrh  |  w)J  =  0. 

(82) 


FIRST  MOMENT  EQUATION: 

Now  we  can  write  the  first  moment  of  (72)  by  left  multiplying  by  {uj\p {  as 
follows 

e2e2  -  - 

e2(w  |  pi  |  dtut)  +  e(u>  |  picCj  |  dju)  +  I  PiLjCk  |  (fdkw)  =  0,  (83) 

where  the  R.H.S.  vanishes  because  of  the  second  constraint  (74)  on  the  collision 
operator.  Because  p  and  £  commute,  adding  (83)  to  its  adjoint  equation  gives 
us  the  following  first  moment  equation 

A  ^2^2  r  A  A  A  A  \ 

e2dt(oj  |  pi  |  w)+edj(u>  |  p^Cj  |  w)+-^-  |  Pi£j£k  \  9jdku>)  +  {dkdju>  |  £k£jpi  |  u>)j  =  0. 

(84) 

The  zeroth  and  first  momentum  equations  (82)  and  (84)  are  partial  differential 
equations  in  the  matrix  elements.  The  macroscopic  equations  of  motion,  a  mass 
continuity  equation  and  a  Navier-Stokes  equation,  come  from  (82)  and  (84), 
respectively. 

We  can  now  determine  the  partial  differential  equations  that  describe  the 
dynamics  of  a  quantum  lattice  gas  in  local  equilibrium.  Inserting  (78)  into  the 
zeroth  moment  equation  (82)  and  retaining  terms  up  to  second  order  in  the 
smallness  gives 

e2dt  |  u/°))  j  +  edi  |  £*  |  =  (85) 

-  ?^-di  (V°)  I  iij-'ij  I  V(0))  +  (0;w(o)  I  \  «<°>)) 

-  (<u>(0)  i  CiCj  |  didjJV)  +  {dtdjjv  |  tit i  I  w<°))) , 
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using  the  fact  that  m  =  ml.  Identifying  the  local  equilibrium  mass  density  (46), 
p  =  |  ««»),  and  momentum  density  (47),  pvi  =  mc(uj(-°'>  |  4  |  u/°l),  we 

arrive  at  the  hydrodynamic  equation  for  mass  flow 

e2dtp  +  edi  (pvt)  -  (86) 

_  I2f^ft  ((jo)  |  IJ-ijr.  |  d.Jo)}  +  (d.jo)  |  £.(j-i)*c.  |  w(o)) j 

-  (V0)  |  44  |  ftftu/0)>  +  (ft  v<0)  I  44  I  «(0)>)  , 

Similarly,  by  inserting  (78)  into  the  first  moment  equation  (84)  and  retaining 
terms  up  to  second  order  in  the  smallness  gives  the  hydrodynamic  equation  for 
momentum  flow 

(pvi)  +  eftng^=  (87) 

-  e2mc2£dj  (<w<°>  |  44  j_14fc  I  ftw(0)>  +  (ft:^(0)  |  4(i~1)*44  I  ^(0)>) 

2 

-  ~m<?£  (<w<°>  1 444  I  4ft^(0)>  +  1 44-4 1  ^(0)>) , 

where  the  ideal  part  of  the  momentum  flux  density  tensor  is  defined  as 

ng“'  =  |  CiCj  |  w<°>).  (88) 

To  obtain  the  macroscopic  equations  of  motion,  we  have  to  determine  the  value 
of  the  matrix  elements  appearing  in  (86),  (87)  and  (88).  This  is  carried  out  in 
Appendix  B.  The  result  is  that  in  the  incompressible  limit  (V  •  v  =  0),  all  the 
matrix  elements  on  the  R.H.S.  of  (86)  vanish  leaving  us  with  a  mass  continuity 
equation 

dtp  +  di(pvi)  =  0,  (89) 

and  the  matrix  elements  on  the  R.H.S.  of  (87)  do  not  vanish  leaving  us  with  a 
viscous  Navier-Stokes  equation12 

dt(pvi)  +  dg(gpviVj)  =  -diP  +  pvd2Vi.  (90) 

In  (90),  the  pressure  is 


and  the  kinematic  shear  viscosity  is 

(92) 

The  form  of  the  kinematic  viscosity  (92)  is  identical  to  that  for  the  classical 
lattice  gas  as  found  by  Henon  [20],  However,  in  the  case  of  the  quantum  lattice 
gas,  the  value  of  Kn  can  be  different  than  the  classical  value  because  of  quantum 
mechanical  interference  of  outgoing  collision  possibilities. 

12As  is  the  case  for  the  classical  lattice  gas,  there  is  a  density  dependent  prefactor  appearing 
in  the  convective  term  and  the  pressure,  g(d)  =  T^T  • 


20 


11  Summary 

In  this  paper,  we  predicted  the  macroscopic  behavior  of  a  lattice-gas  quantum 
computer.  The  following  observations  are  made: 

1.  The  >1' -space  unitary  collision  matrix,  C,  is  successively  blocked,  first  over 
the  on-site  2s-dimensional  manifold,  C  =  0^=1  &  ■  Each  block  U  is  also 
block  diagonal  over  the  equivalence  classes. 

2.  The  projection  operator,  f ,  periodically  causes  a  partial  collapse  of  the 
on-site  superpositions.  Therefore,  a  quantum  computer  with  many  qubits 
can  simulate  a  quantum  lattice  gas  with  only  short-term  and  short-range 
entanglement  and  coherence  of  qubits. 

3.  Streaming  of  complex  amplitudes  (or  the  associated  qubits)  occurs  in  anal¬ 
ogy  to  streaming  in  a  classical  lattice  gas  and  does  not  cause  global  en¬ 
tanglement  in  the  quantum  lattice  gas  because  of  the  application  of  T. 

4.  The  quantum  lattice  gas  can  be  understood  as  existing  between  two  limits, 

a  fully  coherent  quantum  system  and  a  classical  system.  And  its  single¬ 
particle  distribution  function,  fa  =  \uia\2  =  =  (ift\na\ip),  has  the 

form  fa  =  1/  (exp(etp  +  (3ea  ■  p  +  *)E)  +  1). 

5.  Like  the  lattice  Boltzmann  equation  approach  to  simulate  fluid  dynamics, 
the  quantum  lattice  gas  is  a  noiseless  method  that  directly  codes  the  par¬ 
ticle  dynamics  at  the  mesoscopic  scale.  However,  unlike  the  lattice  Boltz¬ 
mann  BGK  collision  operator,  the  quantum  lattice-gas  collision  operator 
obeys  detailed  balance.  Hence,  the  method  is  unconditionally  stable. 

6.  The  macroscopic  hydrodynamic  behavior  is  described  by  a  viscous  Navier- 
Stokes  equation. 

12  Closing  Remarks 

To  mimic  the  behavior  of  other  physical  systems,  quantum  lattice  gases  need 
many  qubits.  A  first  generation  quantum  computer,  with  only  two  qubits, 
cannot  test  the  behavior  of  the  quantum  lattice  gas  at  the  macroscopic  scale. 
However,  useful  tests  could  be  conducted  on  a  network  of  these  first  generation 
machines  to  test  the  practicality  of  the  quantum  lattice-gas  formalism.  For 
example,  we  could  test  the  reliability  and  computational  speed  of  a  network 
of  quantum  computers.  It  is  reasonable  to  expect  that  the  number  of  qubits 
will  grow  exponentially  according  to  Moore’s  law  as  various  quantum  computer 
designs  are  realized  over  time. 

We  know  from  experience  with  classical  lattice  gases  that  although  the  un¬ 
derlying  microscopic  dynamics  is  reversible,  dissipative  shear  viscosity  arises  at 
the  macroscopic  scale — entropy  increases  while  at  the  same  time  information  is 
conserved  because  of  microscopic  reversibility.  The  reason  for  this  is  that  in¬ 
formation,  initially  stored  in  the  spatial  correlations  of  arrangement  of  particle 
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occupancies,  in  time  is  transferred  into  high  order  particle-particle  correlations. 
The  same  process  should  occur  in  a  quantum  lattice  gas.  We  do  not  yet  know, 
in  a  quantum  lattice-gas  setting,  if  there  is  a  way  to  block  this  informational 
transfer  mechanism  and  thereby  reduce  dissipation  at  the  macroscopic  scale. 

Why  consider  constructing  a  quantum  computer  following  the  lattice  gas 
paradigm,  when  a  general  purpose  quantum  computer  could  simulate  a  quantum 
lattice  gas?  The  answer  to  this  question  is  three-fold:  (1)  because  any  member  of 
the  appropriate  unitary  group  associated  with  an  equivalence  class  block  of  the 
collision  operator  is  sufficient  for  the  recover  of  Navier-Stokes  hydrodynamics 
at  the  macroscopic  scale,  so  the  lattice-gas  quantum  computer  is  robust;  (2) 
short-term  coherence  among  only  a  small  number  of  nearby  qubits  is  needed; 
and  (3)  its  behavior  is  predictable  by  analytic  means. 

The  quantum  lattice  gas  method  can  be  straightforwardly  applied  to  three- 
dimensional  fluid  simulations  (the  two-dimensional  case  was  treated  in  this  pa¬ 
per  because  of  its  simplicity)  and  also  applied  to  model  other  physical  systems. 
Lattice  gases  are  a  special  case  of  cellular  automata  where  conservations  and 
detailed  balance  are  imposed  and  an  isotropic  spatial  lattice  is  used.  With 
these  few  restrictions  removed,  the  method  presented  in  this  paper  represents  a 
general  computational  system  called  a  quantum  cellular  automaton. 
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A  Single-Particle  Distribution  Function 

The  single-particle  distribution  function  has  the  form 

(93) 

where  the  natural  log  of  the  fugacity 

In  za=  ap  +  /3ea  -p  +  ~/E  (94) 

is  a  linear  combination  of  the  conserved  scalar  quantities,  the  mass  p,the  mo¬ 
mentum  component  ea  ■  p  along  the  lattice  direction  ea,  andthe  energy  E  at 
a  lattice  site.  The  real  numbered  coefficients  a,  (3 ,  and  7  are  free  parameters 
that  we  will  determine.  It  is  convenient  to  define  the  momentum  and  energy 
independent  part  of  the  fugacity  as 


-Zo 


=  eop. 


(95) 


Since  fa(z0)  =  d  is  the  reduced  density,  d  =  ,  we  must  set 


Zq 


1  -d 
d  ' 


(96) 


This  fixes  the  coefficient  a.  To  fix  the  coefficients  /?  and  7,  we  can  specify  two 
moments  of  the  single-particle  distribution  function  as  constraint  conditions. 
We  begin  by  Taylor  expanding  the  single-particle  distribution  function  f(za) 
about  z0 

f(*a)  =  d  +  f'(za)Sz  +  ±f"(zo)(6z 2)  +  •  •  • .  (97) 

The  derivative  of  /  evaluated  at  za  are 

f(z)  =  (7^2  -»  f(z o)  =  -d2  (98) 

and 

r(z)  =  (7TTF  f"[Zo)  =  2rf3’  {99) 

SO 

f(za)  =  d[l-d6z  +  d?(6z)2].  (100) 

To  determine  5z,  we  begin  by  writing  the  fugacity  in  series  form 


E 

.k= 0 


(0ea-p)k 

k\ 


E 

Lfc=o 


(7  E)k 
k! 


Za  =  Z0 


(101) 
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In  the  subsonic  limit,  p  me,  keeping  terms  only  to  second  order  in  the  velocity, 
the  fugacity  becomes 

za  =  za  1  +  pea  ■  p  +  ^(/3ea  ■  p)2  (1+7 E)  +  0(v3).  (102) 

since  p  ~  v  and  E  ~  v2.  Then  to  second  order  in  the  velocity,  the  change  in  za 
is 

5za  =  za  -  z0  =  flea-p+ ^(flea -p)2 +  yE  +  0(v3)  (103) 

and  the  square  of  the  change  is 

(SZaf  =  (^)V(^a  • P)2  +  0(V 3).  (104) 

Inserting  the  expressions  for  Sz  and  (6z)2  into  the  Taylor  expansion  of  f(za)  we 
have 

f(za)  =  djl-(l-d)  f3ea.p+±(f3ea-p)2+'yE  +  (1  -  d)2  (ea  •  p)2  J 

=  d  [l  -  (1  -  d)  (pea  ■  p  +  7JB)  +  |(1  -  d)(l  -  2 d)p2  (e0  •  p)2 j  .(105) 

We  have  the  freedom  to  choose  the  coefficients  (3  and  7  to  parameterized  the 
distribution  function  as  we  see  fit  to  satisfy  any  two  constraints.  Consider  a 
parameterization  that  fixes  the  value  of  the  coefficients  /?  and  7  by  using  the 
following  moments  for  the  mass  density  and  momentum  density 

B 

m^2fa  (106) 

a=l 
B 

mc^^eafa.  (107) 

a=l 

The  parameterization  may  be  termed  the  non-Galilean  parametrization.  Con¬ 
straints  (106)  and  (107)  are  typically  used  in  the  formulation  of  classical  lattice 
gases.  The  single  particle  distribution  function  using  this  non-Galilean  param¬ 
eterization  was  first  found  in  the  mid  1980’s  by  the  US  researchers  Wolfram 
and  Hasslacher  and  by  the  French  researchers  FYisch,  d’Humieres,  Lallemand, 
Pomeau,  and  Rivet  [1,  21].  Their  derivation  of  (112)  given  below  is  different 
then  the  derivation  presented  in  this  section;  they  used  only  two  free  coefficients 
in  the  expression  for  the  fugacity,  one  for  the  mass  and  the  other  for  the  momen¬ 
tum;  whereas  we  use  three  free  coefficients.  The  reason  for  using  only  two  free 
parameters  is  that  in  the  standard  single-speed  classical  lattice-gas  construction, 
the  energy  is  degenerate  with  the  mass,  so  it  was  deemed  unnecessary  to  keep 


P  = 

pv  = 
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a  separate  free  coefficient  for  the  energy.  Using  (106)  and  (107)  as  constraint 
equations  gives  us  a  non-unity  density-dependent  prefactor  in  the  convective 
term  in  the  hydrodynamic  flow  equation. 

Inserting  (105)  into  (107),  the  odd  term  in  the  distribution  function  expan¬ 
sion  survives  the  first  moment  sum  over  lattice  directions;  the  odd  term  is  the 
one  linear  in  the  momentum.  This  fixes  the  value  of  /?  to  be 

»~~^2  <™> 


so  the  distribution  function  becomes 

fa=d 


n2  j  _  o fi 

1  +  Dea  ■  p +  —  - - —  (ett  ■  p)2  +  (1  -  d)iE 


2  1-d 


(109) 


Inserting  (109)  into  (106),  all  the  even  terms  that  survive  the  sum  over  lattice 
directions  must  add  to  zero.  This  fixes  the  value  of  7  as  follows 

D  1  — 2d 


2  1  -d 


p  —  ( 1  —  d)yE  =  0 


or 


7E  =  D 


1-2  d  p2 


(1  -  d)2  2  ' 

Therefore,  the  non-Galilean  parameterized  distribution  function  is 


fa^d 


1  +  Deaipi  +  -~D  +  2\(d)QnijPiPj 


(110) 

(111) 

(112) 


where  the  density  dependent  prefactor  g(d)  is  defined 

,  ,,  _  D  1-2 d 
9 ^  ~  D+  2  1  -  d 

and  the  traceless  second-rank  tensor  Qa  is  defined 


n  -  =  e  e  •  - 

'vaij  —  catca;  • 


(113) 


(114) 


Qa  is  an  isotropic  symmetric  tensor.  This  mass-energy  degeneracy  leads  to  an 
anomalous  description  of  the  lattice-gas  fluid’s  behavior.  Let  us  see  why.  The 
second  moment  of  (112)  gives  the  momentum  flux  density 


B 

me2  ^2,  €ai&aj  fa  =  P5ij  +  gpviVj.  (115) 

0=1 


The  density-dependent  prefactor  g  appears  in  the  nonlinear  convective  term, 
so  this  parametrization  does  indeed  give  rise  to  non-Galilean  fluid  flow.  The 
pressure  in  (115)  has  a  spurious  quadratic  velocity  dependence 

P  =  pc2  (l  -  g^.  (116) 
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k  » 


B  Determination  of  the  Matrix  Elements 


For  the  triangular  lattice,  tensors  made  up  of  products  of  the  lattice  vectors  are 
symmetric  and  isotropic  [1],  We  have  the  following  identities 


(1  |  Ci  1  1)  =  0 

(117) 

£-iCj  \  l)  =  jjSij 

(118) 

|  CiCjCk  I  1)  =  0 

(119) 

D(D  +  2)  ^'^kl  +  ^ufyk) 

(120) 

(1  |  CiCjCkCt  |  1) 

The  Jacobion  of  the  collision  operator  is  circulant.  Its  eigenvectors  correspond¬ 
ing  to  the  nonzero  eigenvalues  span  the  kinetic  space ,  which  contains  a  viscous 
subspace  characterized  by  a  degenerate  eigenvalue,  denoted  by  kv  [22].  The 
eigenvectors  in  the  viscous  subspace  are  j CiCj  |  1),  for  i  ^  j.  Therefore,  we  have 

JCiCj  I  1)  =  KrjCiCj  I  1),  (121) 

or  inverting  this  over  the  kinetic  viscous  modes 


J-'CiCj  |  1)  =  — CiCj  1 1). 


(122) 


Using  these  identities  along  with  the  epsilon  expansion  of  |  u>),  we  can  work 
out  the  value  of  the  matrix  elements  which  appear  in  the  mass  and  momentum 
hydrodynamic  equations  of  the  quantum  lattice  gas  at  the  macroscopic  scale. 
We  have 


<o/°>  |  CiCj  |  u>(°>) 


d(  1 1  CiCj  1 1)  +  ^  (1  |  CiCiCu  |  l)t>*  + 

gdD(D  +  2) 

2c2 

dB x  t  gdB  lx  r  ,  r  r  ,  f  r  ..  gd(D  +  2)Bs.  2 

j-)  *ij  '  2^2  '  OifcOjl  ~r  OilOjk)  2,D(?  VijV 

„2\  _ 

(123) 


(1  |  CiCj  (ckCt  -  5j^j  |  l)vkVl 


dB  (  v2\  „ 


(0) 


<«(0)  I  CiCjCk  |  djdku>W) 


CiCj  I  didjU>W)  =  ^(1  I  CiCjCk  I  l)d,djVk 
=  0. 

dD 


(124) 


2c 


(1  |  CiCjCkCi  |  1  )djdkvi 


dB 


2  c(D  +  2) 
dB 

2  c(D  +  2) 


(5ij6ki  -(-  5ikSj,  4-  SuSjk)  djdm 


(2didkvk  +  d2Vi)  . 


(125) 
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<u,(0)  I  Cij-'Ci  I  vm)  =  2^(1 1  CiJ-'Cjtk  I  i)a,-«t 

=  1^(1 1  C-iCjCk  1 l)c^ 

=  0.  (126) 

(w(0)  I  I  dkL>W)  =  ^(1 1  CiCjj-'CkCt  I  i)afcVl 

=  2e^  (1  |  CiCjCkCi  I  l)dfcv( 

=  2^ZTT2}  (^'4'  +  ^  +  5il6jk)  dkVl 
=  2^T2)  ( +  ^  +  «  •  (127) 
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